A Chain-Boson Model for the Decoherence and Relaxation 
of a Few Coupled SQUIDs in a Phonon Bath 
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We develop a "chain-boson model" master equation, within the Born-Markov approximation, for 
a few superconducting quantum interference devices (SQUIDs) coupled into a chain and exchang- 
ing their angular momenta with a low temperature phonon bath. Our master equation has four 
generators; we concentrate on the damping and diffusion and use them to study the relaxation 
and decoherence of a Heisenberg SQUID chain whose spectrum exhibits critical point energy-level 
crossings, entangled states, and pairs of resonant transitions. We note that at an energy-level cross- 
ing the relevant bath wavelengths are so long that even well-spaced large SQUIDs can partially 
exhibit collective coupling to the bath, dramatically reducing certain relaxation and decoherence 
rates. Also, transitions into entangled states can occur even in the case of an independent coupling 
of each SQUID to the bath. Finally, the pairs of resonant transitions can cause decaying oscillations 
to emerge in a lower energy subspace. 

PACS numbers: 03.65. Yz, 03.67.Mn, 03.67.Pp, 75.45.+j, 74.50.+r, 85.25.Dq 
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I. INTRODUCTION 

A. SQUIDs and SQUID Chains 

A superconducting quantum interference device 
(SQUID) can be made from a small strip of aluminum 
bent into a ring, joined at the ends, and cooled to a 
milliKelvin temperature. Aluminum is a superconductor 
and an aluminum oxide layer, where the ends meet, forms 
a "Josephson Junction" potential barrier. Precise tuning 
of an externally imposed magnetic flux can cause oscil- 
lations of current^ between clockwise "|t)" and coun- 
terclockwise "|D" states via an evolving phase inter- 
ference, e i( E i-Eo)t/h ^ b e t ween the ground state, |0) = 
(lT) + U»/\/2, and first excited state, |1) = (|T)-|-l))/V2, 
which span the low-energy dynamics of the device and 
form a logical basis for a "qubit" of quantum informa- 
tion, Cq |0) + Ci 

Many SQUIDs can be coupled together into a chain. 
The aluminum rings are not topologically linked but their 
proximity allows capacitive and inductive interaction be- 
tween nearest neighbor SQUIDs-^. We are principally 
interested in using the chain to encode 5 -^ 7 - and protect 
or correc t 8 : 9 : 10 ' 11 quantum information. "Logical gate" 
operations on a chain of SQUIDs can create^ 2 - useful 
non-classical correlations quantified by their "entangle- 
ment of formation"—. Chains can propagate excitations, 
qubit a 14 ' 15 : 16 ' 17 , and even entangled singlets, (|t|) — 
They could algQ provide lon 

g-sought ex- 
perimental realizations of spin chains which in turn illus- 
trate the correlations and phase transitions of many-body 
physics. Cooling a chain to its zero-temperature ground 
state can prepare useful entanglement 21,22 . Ground state 
entanglement can vanish abruptly, for example, as spin- 
spin couplings are adjusted across a critical poin t 23 : 24 ' 25 . 
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FIG. 1: Relaxing through entangled states. A Heisenberg 
two-SQUID chain above the critical point B c cooling from 
1 11) to 1 00} . In the process, the states 1 01) ± 1 10) are oc- 
cupied, resulting in a surge of the entanglement of forma- 
tion, even though the ground state is separable and the large 
SQUIDs are dissipating independently. (8J/h — 1.0GHz, 
2B/h = 1.5GHz, k B T/h = 0.3GHz, R = 10/im, I = 3^A, 
p — 5g/cm 3 , c± — 5km/s). 



In that case, it is the intermediate energy states that can 
possess entanglement; their quantum correlations can be 
evident, at a warmer but not too-high temperature, when 
mixed sufficiently into the equilibrium stat o 26 ' 27 . 



B. Decoherence and Relaxation 

Here we consider some effects of exposing a SQUID 
chain to a low temperature bath of phonons. We suppose 
the SQUIDs are lithographically etched and deposited 
into a solid crystal (e.g. of silicon) which, for simplicity, 
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we assume surrounds the chain. In a SQUID 's oscillation 
between current states, the conservation of total angular 
momentum requires torsional oscillations of the solid and 
thus the emission and absorption of phonons 2 ^. Deco- 
herence of superconducting qubits is also caused by two- 
level defects in the Josephson Junction barrier material 29 , 
which could in principle be purified, whereas here we con- 
sider only the minimal unavoidable coupling required by 
symmetry. 

Usually, when a quantum system is opened-up to its 
finite-temperature environment, its energies and eigen- 
states are perturbed and the new energies and station- 
ary states are viewed as renormalized quantities. Then it 
equilibrates, by decoherence and relaxation, to a station- 
ary thermal mixture of these eigenstates. In the relax- 
ation process, the populations (i.e. probabilities) of each 
eigenstate are adjusted until a thermal mixture is ob- 
tained. Decoherence is the decay of the oscillations (e.g. 
of the current) . Some decoherence goes along with the re- 
laxation — adjusting the probabilities of the eigenstates 
undermines the support for any phase interference be- 
tween them — but it can also be caused by pure dephas- 
ing, in which the phase's probability distribution spreads, 
reducing the average phase interference, or "coherence," 
without adjusting the populations. As described in ap- 
pendix [3 the populations and coherences are, respec- 
tively, the diagonal and off-diagonal matrix elements, in 
the energy eigenbasis, of the chain's density operator pit). 
The absolute value of a coherence is just the (decaying) 
envelope of the oscillations. 

We are interested in large SQUIDs of radius R — 10/im 
with a 3[iA current oscillating at 1.0 ± 0.5 GHz. For a 
single SQUID the decoherence model gives decoherence 
times of a few /xs2^ which are not inconsistent with those 
of recent experiments^^. On the other hand we will also 
simulate small SQUIDs, with R = 10 nm and I = 0.1//A, 
to demonstrate the concept of a collective coupling to 
the environment. The concept of a collective coupling 
is relevant to the critical point open system dynamics of 
the large SQUIDs. 



C. The Chain-Boson Model 

The spin-boson model 3 - 1 - has been widely applied to- 
wards a better understanding of the environment's effect 
on a single qubit's coherence. In the spin-boson model, 
a two-level system is coupled to an environment of os- 
cillators which model a heat bath, such as is often used 
in studies of quantum Brownian motio n 32 i 33 ' 34 . Because 
of the coupling, the system becomes entangled with the 
bath. When averaged over environmental outcomes, the 
system typically loses coherence and thermalizes. Chains 
of harmonic oscillators interacting with a finite temper- 
ature bath have been considered m 35 i 36 i 37 . 

In light of qubit encoding schemes and the desire to 
process and protect quantum information, it is necessary 
to study the decoherence of multiple qubits. In the chain- 
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FIG. 2: Decaying oscillations emerging in a lower energy 
subspace. Three large SQUIDs dissipating independently. 
\ip4.) + 1^7) cooling through a mixture including + ^4). 
The transitions |t/>4) <-» |t/)i) and \ij)r) «-> IV^) are resonant. 
P47, which decays ~ e -1 " 47 ', is absorbed into P14 which decays 
more slowly ~ e~ Fl4t . It might be easier to start with |TTt), 
close to the equal superposition of |i/>4) = 1 001) + 1 100) + |010) 
and IV'7) = |011) + 1 1 10) + 1 100) , to prepare the coherence 
p 4 7 which flows into p 14 . (6J/h = 1.0GHz, 2B/h = 1.5GHz, 
k B T/h = 0.3GHz, R = 10^m, I = 3/iA). 

boson model one embeds a chain of qubits in a bosonic 
bath so that the qubits experience a location- dependent 
interaction with the bath variables. 

For a system comprising a register of qubits, two types 
of system-bath coupling have already been extensively 
considered. The simplest is a collective coupling, in which 
each qubit couples to the same environmental variable 
as the rest. This is appropriate to scenarios where the 
qubits are spaced closer together than the relevant wave- 
lengths of the bath, i.e. those corresponding to the qubits' 
transition frequencies. The symmetries of the collective 
coupling can lead, with certain system Hamiltonians, to 
decoherence free subspace o 10 ' 11 ! 38 and Dicke super- and 
sub-radiance and super- and sub-decoherenc o 39 ' 40 . 

The other commonly used type of system-bath cou- 
pling is the independent coupling model, in which each 
qubit couples to its own bath, separate from the baths 
used for the other qubits. This is an appropriate model 
for qubits spaced farther apart than the relevant wave- 
lengths of the bath. And in the context of solid state 
qubits, a significant source of noise is the voltage leads 
that control the qubits. With one lead per qubit, the 
independent baths model is a natural assumption. The 
independent coupling model has been used to study the 
decoherence during two-qubit logic gate o 41 ' 42 as well as 
the entanglement rate for coupled qubit o 22 i 43 ' 44 . 

When there is one voltage lead controlling multi- 
ple qubits, one typically uses the collective coupling 
model, as each qubit is experiencing the same electronic 
noise. A likely scenario for a pair of qubits is an in- 
dependent lead for each qubit as well as one common 
lead. The disentanglement and decohering of the two 
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qubits in two separate and one common cavity have 
been considere d 45 ' 46 ' 47 ' 48 ' 49 ' 50 . Several works compare 
the collective and independent bath scenarios for cou- 
pled qubit a 51 ' 52 ' 53 ' 54 . Fine-tuning the inter-qubit cou- 
pling to protect against collective dissipation has also 
been studie d 55 ' 56 . In comparison, the optimum qubit- 
qubit coupling was examined for the case of independent 
dissipatio n 57 ' 58 . 

Most studies fall into these two categories, whether 
the qubits are coupled together or not. But for uncou- 
pled qubits there has been careful consideration of the 
intermediate scenario, in which the qubits are neither 
far apart nor close togethe r 59 ' 60 ' 61 . The relaxation and 
decoherence rates depend on k ■ fjk = wTjk, where the 
k are the bath wavevectors that interact with the reg- 
ister of qubits at its transition frequency cj, and Tjk is 
the phonon transit time between qubits j and k sepa- 
rated by fjk- These results make an elegant transition 
between the two limiting cases of qubits close together 
and far apart. In the intermediate scenario the bath can 
induce entanglement between uncoupled qubits^ 2 -, as can 
also happen with a collective coupling to the bath—. Un- 
coupled SQUIDs can also be entangled when the bath is 
specialized to a single mode cavity^ 4 -. 



D. Role of the Critical-Point 

In our chain of SQUIDs coupling to phonons, an es- 
sential point is that the inter-qubit couplings can play a 
crucial role in determining which are the relevant bath 
wavevectors and thus whether the qubits couple collec- 
tively or independently to the bath. In our model, there 
is a critical point energy-level crossing at which some 
transition frequencies are so slow that k ■ fjk — u) Tjk <C 1 
even though the qubits may be far apart with respect 
to the uncoupled qubits' transition frequencies. At the 
critical point the chain can thus obtain some or all of 
the benefits of a collective coupling to the bath. Some 
interesting studies include the intermediate regime for 
coupled qubit a 65 ' 66 and detail the sub-radiant behavior 
of dipole-coupled qubits for small k ■ fj^, but the role of 
the inter-qubit coupling in determining the relevant bath 
wavevectors and their relation to the inter-qubit distance 
is not emphasized. 



E. Overview 

Working from the total Hamiltonian H = Hs+Hb+V, 
describing the SQUIDS, the bath, and their coupling, 
we apply master equation techniques from the quantum 
Brownian motion mode l 32 ' 33 ' 34 , in the Born-Markov ap- 
proximation, to a chain of SQUIDs interacting with a 
phonon heat bath via location-dependent couplings. The 
resulting generators of the open system dynamics are as- 
sociated with four types of coefficients: the renormaliza- 
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FIG. 3: Decoherence and Relaxation in two large SQUIDs dis- 
sipating independently. 1 00) + |11) cooling above the critical 
point, (a) The population dynamics are the same for an equal 
mixture of \ipi) = |11) and |?/>4) = 1 00) ; only the superposition 
is entangled, (b) The entanglement does not oscillate because 
any phase between 1 00) and |11) may be generated locally. It 
decays faster than the coherence pi4 ~ e~ Fl4t upon which it 
depends. Later, the populations mix back in some entangle- 
ment. (8J/h = 1.0GHz, 2B/h = 1.5GHz, k B T/h = 0.3GHz, 
R = 10pm, / = 3fiA). 



tion, anomalous diffusion, damping, and diffusion. We 
use the damping and diffusion to develop a matrix ele- 
ment equation for the populations and coherences. It is 
similar to a Bloch-Rcdficld 68 equation but in the inter- 
action picture. It gives the relaxation and decoherence 
rates, as well as the possibility for coherent oscillations to 
moves 9 - from one subspace, where they are decaying, into 
a lower energy subspace with a longer decoherence time. 
A complete network of selection-ruled transitions leads 
to thermalization. This is usually the case for large well- 
spaced SQUIDs. For small closely-spaced SQUIDs the 
network is broken by the degeneracies of their collective 
coupling to the bath; decoherence and transition rates 
can scale with the number of SQUIDs or vanish; some 
subspaces are protected. Finally, at the critical point 
even the large well-spaced SQUIDs can acquire some of 
this collective Dicke super- and sub-radiant and super- 
and sub-decoherent behavior. 



II. THE CHAIN, BATH, AND THEIR 
COUPLING 



A. The Heisenberg SQUID Chain 

For its mathematical simplicity and relevance to quan- 
tum information processing, we consider the isotropic 
Heisenberg coupling between nearest neighbor SQUIDs. 
In principle this can be engineered with a precise balance 
of inductive and capacitive coupling between nearest- 
neighbor SQUIDs. However, several of our methods are 
applicable to other types of coupling. The "antiferromag- 
netic" (J > 0) Heisenberg chain of N SQUIDs evolves by 
its Hamiltonian, 



H, 



N 



Bo*), 



(1) 
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assuming periodic boundary conditions <tjv+i = <?i- 
Here, aj are the Pauli matrices for the jth SQUID, with 
|t) and ||) the eigenstates of a z , and B is an energy split- 
ting (not an imposed magnetic field) appropriate to the 
natural precession of our pseudo-spin qubits: coherent 
quantum oscillations between ||) and ||). 

The Heisenberg and magnetic sums commute and their 
respective quantum numbers I and m determine the en- 
ergy spectrum {IJ — mB} (up to some degeneracies not 
split by J and B). Regardless of N, each eigenstate is a 
linear combination of states with the same number of |1) 
vs. |0) qubits (m = Nq — Ni) and is typically entangled, 
with the exception of the extremal m states |11 . . .) and 
|00 . . .). Increasing B relative to fixed J causes energy- 
level crossings. At a critical point, B c , the ground state 
changes from entangled to unentangled. 

The two-SQUID chain's energies are shown in Fig- 
ure Q] as a function of the single-SQUID energy splitting 
B. Figure [2] plots the energies for a three-SQUID chain. 
Their respective eigenstates are as follows: 




FIG. 4: Entanglement oscillations in two large SQUIDs dis- 
sipating independently. |01) cooling below the critical point, 
(a) The population dynamics are the same for an equal mix- 
ture of |V> 2 > = |01> 1 10) and|^ 3 ) = |01) + |10). (b) From their 
initial superposition |01), the evolving phase between \tp2) and 
\tp3) drives « 6 x 10 5 rapid oscillations of the entanglement of 
formation (eof); we have plotted a moving average <eof>, the 
upper bound (>eof), and the lower bound (<eof) which is the 
same as the entanglement of the mixture. Initially, the mov- 
ing average decays with the coherence P23 ~ e -1 " 23 * but later 
pulls away to equilibrate. (8J/h = 1.0GHz, 2B/h = 0.5GHz, 
k B T/h = 0.1GHz, R = 10/ra, / = 3/j.A). 
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B. The Harmonic Phonon Bath 

The harmonic crystal Hamiltonian is composed of 
phonon modes labelled by wavevector k and polarization 
index s. Each phonon contributes an energy huj s (k): 

H B = J2hu; s (k)ai s a~ ks . (2) 

fc,s 

Here at and a? are the phonon creation and annihila- 

ks ks * 

tion operators. For frequencies below a "Debye" or cut- 
off frequency A we assume linear dispersions, uj\^{k) = 
c± \k\ and u>3(k) = cm |A;|, for transverse (_L) and longitu- 
dinal (||) polarizations e s (k) = |fc±i, k±_2, k\\ \- We will 
also make use of the momentum density at site r: 

= -i g %u s {k)p Hs& ~^* e e s (k) (3) 

ks 

with V and p the volume and mass density of the crystal. 



C. The Chain-Bath Coupling 

The coupling between the SQUIDs and the crystal 
arises from the fact that each SQUID 's current la? is 
formed from the electronic band states in the reference 
frame co-moving with the lattice sites during the tor- 
sional oscillations of the crystal—. In the lab frame the 
electron velocity v e must include the speed u = ir/p of 
the lattice sites: v e — j/(en e ) + u. Here j is the cur- 
rent density and e and n e are the electron charge and 
number density of the electrons. Their kinetic energy 
density n e ^m e \v e \ 2 thus acquires a cross term and the 
total Hamiltonian must include an additional 

V = (m e /e) J d 3 rj-u. (4) 

In appendix [Bl we first derive, following 28 , the coupling 
Vj of an individual SQUID to the crystal and then, be- 
cause the current density j is the sum of the individual 
densities jj, sum their contributions at SQUID locations 
Xj (spaced a distance d apart along the rc-axis) to obtain 

V = \Y,Xj<t>j (5) 
3 

akin to a quantum Brownian oscillator's bilinear cou- 
pling to an oscillator bath, Xx<j), (only here we 
are summing over several contact-points between the 
chain and the bath). The coupling constant A = 
27r/(m e /e) yj Rhc± /2pV fixes the strength of the interac- 
tion between chain operators Xj = crj and bath operators 

4> 3 = ^\RJi(k x R) (afce^ + ay-*' 3 *) , (6) 

k 

for which we have defined the first transverse polarization 
ei (k) to lie in the plane of the SQUID ring and Ji (k x R) = 
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Ji(\k\Rsm6) is the first order Bessel function (and is 
the polar angle between the z-axis of the ring and the 
wavevector k). 

The selection rules, (l'm'\ Xj \lm) ~ 6 m ',m±2 reflect the 
fact that the "interaction operator" Xj simply flips the 
j th |0) or |1), giving a non-zero probability only to obtain 
an eigenstate with one more |1) or |0). The allowed bath- 
driven transition frequencies are u = [(I — I') J ± 2B]/h, 
and transitions driven between degenerate states are thus 
between distinct I and I' at a crossing of energy-levels. 

The two-SQUID chain's interaction operators are, in 
its energy eigenbasis, 



(a) Populations ^ Coherences 




Xl.2 = 



V2 



( Tl 1 \ 
Tl ±1 
10 1 

V o ±11 o / 



(7) 



They present a "network" of selection-ruled transitions: 
\ip2) \ipi) and \tpi) <->■ \ips) «-> \ipi) with transi- 
tion energies 2B and 2B ± 8 J. 

For future reference, this network should be com- 
pared with the case of a collective interaction operator 
J z = J2j Xj in which (for two SQUIDs) the ±1 matrix 
elements above cancel and the selection-rules constrict 
to \ipi) <-» 1^)3) <-> ^4) leaving \ip 2 ) as a protected (1- 
dimensional) subspace. 



FIG. 5: Dicke Super- and Sub-radiance and Super- and Sub- 
decoherence in two small SQUIDs dissipating collectively. 
I^i) + 1^2) + 1^3) + IV'*) cooling above the critical point, 
(a) 1^)4) flows into \tp3) which flows into \ipi). The singlet 
|^2) = 1 01) — 1 10) is a protected subspace and cannot relax 
into the populations do not approach their thermal levels 
(gray-lines). But the collective relaxation rate from \1jj3) into 
\ifii) is double the independent rate, (b) The coherence P12, 
between and \ip2), barely decays; the only mechanism 
for that is the \ipi) — > |i/>3) transition which is suppressed by 
the cold bath. Note also how P34 flows, in a pair of reso- 
nant transitions, into pn whose collective decoherence rate 
is double its independent dissipation rate. (8J/h = 1.0GHz, 
2B/h = 1.5GHz, k B T/h = 0.3GHz, R = Warn, I = O.lpA). 



III. FORMALISM 
A. The Born Approximation 

We now apply master equation techniques from the 
model of quantum Brownian motio n 32 i 33 i 34 to the chain 
of SQUIDs interacting with their phonon bath. 

The unperturbed Hamiltonian Hq = H$ + Hb de- 
fines a standard interaction picture in which the coupling 
V(t) = e iHot/hy^ e -iH t/h determines the evolution of 
the system-bath density operator w(t). To second order 
in the coupling its instantaneous rate of change is 

ib = -l.[v(t)Mo)}-^ f & m), iV(t'),w(0)}}- (8) 

Next, we assume the bath is in a thermal state pb — 
(q-Hb /k B T^ 1 <2 IR initially uncorrelated with the chain: 
w(0) = p(0) (£> pb- Here T is the temperature and 
Z B = J2 n eT En l kBT is the bath's partition function. We 
average (trace) over the states of the bath to obtain a 
preliminary equation of motion for the system: 




(the first order term vanished because Trg[<^>j pb] = 0). 



For a given (j, /c)-pair of SQUIDs, the trace of the 
nested commutators results in four integrals. They differ 
only in the permutation of terms; a representative one is 

rt , * , 

Xj(t) / dt' X k (t')Tr B [^(t)4> k (t')pB]p(0), (10) 
Jo 

V v ' 

where the over- and under-braces highlight the time- 
averaging of the interaction operator X k {if), into what 
we call the noise Vjkif) and susceptibility Ujk{t) oper- 
ators, weighted by the kernels Vj k (r) and p,jk(r) that, 
with t = t — t', are the real and imaginary parts of the 
bath correlator Tr g[(f>j(t) <fik(t') Pb]- (In the context of 
quantum Brownian oscillator systems, the susceptibility 
kernel used here is called dissipation t? 3 . 3 - there.) 

The four integrals for each SQUID-pair, originating 
from the nested commutators and each contributing a 
noise and a susceptibility, can be collected to obtain a 
"Born" (but not yet Born-Markov) master equation (still 
in the interaction picture, still only valid to second order), 

A 2 N 

p-^E (&,\yjk,fl]-wAtiik,p}]) , (11) 

j,k=l 

where we have used Born's approximation that replaces 
p(0) with p(t) w p(0). 
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(a) Populations (b) Coherences 




FIG. 6: Critical behavior of two large SQUIDs. \ipi) + 
l^) + 1^3) + 1^4) cooling at the critical point. The dis- 
tant (e.g. 40/im) SQUIDs are closely-spaced compared to the 
bath wavelengths corresponding to the frequency U12 = 0; 
the <-> \ip2) transitions are blocked, akin to the collec- 
tive behavior of two small SQUIDs close together, (a) \tp4) 
flows into 1^3) and \ip2) = |01) — |10), whereas ^3) flows only 
into giving it an excess population over l^fe) which the 

chain cannot quickly resolve; the \ip2) — * \if>i) — * \i/>s) — » 
pathway takes an extraordinarily long time (not shown), (b) 
The coherence p\2 barely decays; transitions \tp2) <-> \"4>i) are 
blocked while transitions from \ipi) and \ip2) are suppressed 
by the cold bath. Some of the decaying P34 is absorbed by 
pis. (8J/h = 1.0GHz, 2B/h = 1.0GHz, k B T/h = 0.2GHz, 
R = Wpm, / = 3pA). 

Born's approximation adds an updating to the inte- 
grated solution for p(t); in a small time step dt the in- 
stantaneous change dp depends not on the initial p(0) 
but on the updated instantaneous p(t). We will see that 
this updating is needed for the long-time equilibration 
to thermal equilibrium. The Born master equation, al- 
though technically still only valid to second order, is a 
plausible guess at the longer-time open system dynamics 
which, like any other theory, can only be supported by 
real data or exactly solvable open systems, like the quan- 
tum Brownian oscillator, and may not be valid in every 
case. 



B. The Markov Approximation 

Next, we time-average the interaction operators Xk(t') 
into the noise Vjk(t) and susceptibility Ujk(t) operators 
by making a Fourier expansion of the interaction opera- 
tors using amplitudes Xg(t) and V k (t): 

X k (t-T) = J2 {^(t) cos(wr) - V%{t) sin^r)) . (12) 
uj>o 

(In the eigenbasis of Hs, the non-zero matrix elements 
of Xjf are just those matrix elements of Xf. with en- 
ergy difference ±Huj; multiplying those same matrix ele- 



ments by =Fi gives V%. (This sign convention is consis- 
tent with the matrix elements of a harmonic oscillator's 
position, (n\ a) + a \n + 1) = \/n + 1, and momentum, 
i (n\ — a \n + 1) = —i\]n + 1.) X? consists primarily 
of the diagonal matrix elements of Xk but also includes 
matrix elements between degenerate states.) 

We thereby obtain basis-independent expressions for 
the noise, 

V jk (t) = E (D%{t)X£{t) - A%{t)f%{t)) , (13) 

LO>0 

and susceptibility, 

= E - . (14) 

where the coefficients (to be discussed shortly) of diffu- 
sion D, anomalous diffusion A, renormalization r, and 
damping 7 serve to Fourier-sample (at least for t — ► 00) 
the real (vjk) and imaginary (pj k ) parts of the bath cor- 
relator: 

D jk (*) = Jo dTV ik (r) cos ujt Af k (t) = fodTV jk (r) sin ujt 
r jk (*) = Jo dT V 3 k (t) cos ujt 7$. (t) = J^dTp jk (r) sin uit. 

Finally, we make the Markov approximation, that uses 
constant coefficients obtained in the limit t — > 00 in place 
of the time-dependent ones, and substitute these noise- 
and susceptibilty-operators into Eq. (jlip to obtain our 
Born-Markov master equation, 

\2 N 

j,fclu>0 

-D^ k [X 3 , [X^p]]+A^[X 3 , [V^p]]). (16) 

This is essentially a Fourier-series version of the Born 
master equation using Markov (constant) coefficients. It 
is reassuring that, in the limit of only one contact point 
(no sum over j, k) and only one energy splitting (no sum 
over uj) it becomes 

p^ir[x,{x,p}]-iX[x,{p,p}]-D[x, [x,p]]+A[x, [p, p}} (17) 

which is the well-know n 33 ' 34 Born-Markov equation for 
the quantum Brownian motion of an oscillator system 
(Hs — huia^a) with a bilinear coupling (V ~ x<j>) to an 
oscillator bath. 



C. Cross Term Coefficients and Spectral Densities 

In appendices O iDl and [El we study the bath correlator 
and obtain the coefficients of diffusion 

D% = |j, feM coth(^L) = lim D%{t) (18) 
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and damping 

7^ = |^H = t Hm 7 ^(t), (19) 

written in terms of spectral densities Jjk(u>) with, e.g., 

, , hmjl 2 R 4 ui 5 „ .,3, ,5 , . .,. . , 
Jjj(tj) = l^yT^«2 } ' { 2' 3}rM) (20) 

the individual spectral density for a single (j = k) 
SQUID. Here p F q is the generalized hypergeometric func- 
tion and tr = R/cj_ is half the time it takes a phonon to 
traverse the SQUID. 

We then find the same spectral density to first order 
in small k ■ fjk — to Tjk <C 1 for a distinct (J ^ k) pair 
of SQUIDs separated by fjk, where k is a bath wavevec- 
tor at the chain transition frequency uj and Tjk — rjk/c± 
is the phonon transit time between the SQUIDs. Re- 
markably, even "large" and "far-apart" SQUIDs can have 
wTjfc <C 1 at a level crossing where w«0. In these cases, 
as for close-together small SQUIDs, we thus use 

Jjfe(w) = Jjj{u>) when coTjk <C 1. (21) 

On the other hand, the cross term spectral densities 
vanish for truly far-apart SQUIDs satisfying k ■ fjk — 
LJTjk > 1: 

JjkifjJ) = when us Tjk 3> 1. (22) 

We avoid the intermediate regime by tuning the 
chain Hamiltonian close-to or far-from any critical point 
energy-level crossings. 

In our Born-Markov master equation (|16p . the sum 
over transition frequencies thus splits into two parts, one 
for LOTjk 1, for which all cross terms vanish, and one for 
<-<JTjk 1, for which the cross terms are equally weighted. 
For example, the diffusion term becomes 

j.k w>l/r jfc j,k Lj<^l/T jk 

(23) 

The case cuTjk < 1 is the provenance of a collective 
coupling to the bath V = MYlj Xj)4>i (whose square in- 
cludes all cross terms equally). In our case, J2j %j — 
J2j a j = Jz- A collective coupling can have degener- 
ate subspaces with which the system Hamiltonian may 
(or may not) cooperate, giving a decoherence- and/or 
relaxation-free subspace (or not) 11 ' 38 ! 70 . On the other 
hand, when we include all the cross terms equally, the 
master equation has N times as many terms, for N 
SQUIDs, as in the case of independent dissipation. This 
can scale the decoherence and relaxation rates linearly 
with the number of subsystem a 39 i . 




FIG. 7: Automatic quantum error correction in three 
small SQUIDs dissipating collectively. \1jj5} + \^pe) cool- 
ing to \ip2) + \ips}- ( a ) The collective interaction operator 
J z = X\ + X2 + Xz breaks the network of allowed transi- 
tions into three isolated pieces: <-» \1jj2}', \tpa) <-*• l^te); 
and \ips) *-» \ipr) <-> ^4) «-* l^i)- (b) The p56 coherence 
flows along with the populations into P23 because u>25 = ^36- 
Then p~23 never decays; transitions back into ^5) and \ipe) 
are suppressed by the cold bath and are resonant (and there- 
fore coherent) anyway. (6J/h = 1.0GHz, 2B/h = 1.5GHz, 
k B T/h = 0.3GHz, R = lOnm, I = 0.1/LtA). 



IV. THE FOUR GENERATORS 

A. The Renormalization and Anomalous Diffusion 

The master equation has four generators, one for each 
of the coefficients. It can be showr* 7 -^, that for a weak 
coupling of our Heisenberg SQUID chain to the bath and 
with J and B set close-to or far-from any resonant pairs 
of allowed transitions, the renormalization and anoma- 
lous diffusion contribute effectively Hamiltonian dynam- 
ics, which can be dropped from the master equation by 
"renormalizing" the chain Hamiltonian. 

We would like to concentrate on the damping and dif- 
fusion. To make concrete sense of these two generators, 
our examples use the familiar and useful energies and 
eigenstates of the Heisenberg SQUID Chain. In other 
words, we assume that the renormalized Hamiltonian is 
given by Eq. (J!]) as though the "bare" SQUIDs were 
engineered with a slightly different original Hamiltonian, 
including a "counter term" that cancels the renormal- 
ization and anomalous diffusion caused by their crystal 
environment. This engineering may actually be quite dif- 
ficult, relying perhaps on many trials and errors. In any 
case, our general discussion of the damping and diffusion 
is in terms of the renormalized energies and the matrix 
elements in the renormalized eigenbasis. 



B. The Damping and Diffusion 

We now show that the damping and diffusion work to- 
gether to effectively decohere and thermalize the system 
within each network of allowed transitions. In Eq. (|16p 
we set r^ k = = A^ k in accordance with the aforemen- 
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(a) Populations (b) Coherences 




FIG. 8: Critical behavior of three large SQUIDs dissipating 
independently, ttl) cooling at the critical point. The distant 
SQUIDs are closely-spaced compared to the bath wavelengths 
corresponding to the frequencies LU12 = W13 = u)23 = and 
W45 = ^46 = ^56 = 0. (a) The |i/> 2 ) <-» an d l^s) <-» |V>i) 
transitions are cut out of the network while upward tran- 
sitions are suppressed by the cold bath; therefore \ip2) and 
l^) have long-lived populations, (b) P23 is long-lived be- 
cause the transitions out of its ^2) and ^3) supports are so 
rare. (6J/h = 1.0GHz, 2B/h = 1.0GHz, k B T/h = 0.1GHz, 
R = 10/xm, I = 3pA). 



(uj/3a ~ top a ) and states |<5) near E$ (ujpg w tups) into 
those states \j3) that are accessible by transitions from 
I a) and \S) respectively. 

In the first and fourth terms, the remaining transitions 
from states near E a (u)p a ~ up a ) and E$ {uJr$ ~ ojps) 
tend to decrease p a $ in proportion to p a g as well as in pro- 
portion to nearby (in energy) matrix elements p & 8 (in the 
same column) and p a g (in the same row). In the second 
and third terms, the remaining nearly-resonant pairs of 
transitions from \/3) and I7) into \a) and \S) (ui a p w ljsi) 
tend to increase p a s in proportion to those matrix el- 
ements pfj~ t within, and/or nearly within, the diagonal 
that includes p a $ (when the density matrix is stretched to 
be linearly spaced with increasing energy). When S = a 
this diagonal is the central diagonal, sometimes called the 
diagonal. 

C. Decoherence and Relaxation 



tioned renormalization and obtain 



P = 



J,k,u> 



(p» k [Xi , [X£ , p}] +i^ k [Xi , {Vf k ,/>}]). (24) 



We expand out the commutators, work in the energy 
cigenbasis (where Jkv a p = Ep — E a is the energy lost 
in the transition, at frequency cu a p, from \(3) to \a)) by 
inserting resolutions of the identity, e.g. / = J^p \P) 
and define the rates 



r aP 
L jk 



T2 77 JjkiVap) [ 1 + COth( TT7 - 77 ;) 



h 2 2 



•2k B T' 



(25) 



to obtain a matrix element equation, in terms of the pop- 
ulations p aa = (a\p\a) and coherences p a s = ( a \p~\d), 
for their rates of change, 

PaS = - X! ( X! XjapXkPaFjZ XjapX^sT^pp^ 
j,k Pa 0-y 

- ^2 XkapXjfSFjk PPi + ^2 X k5pXjP8~Ff k p a s), (26) 

Pi sp 

where we have used the evenness of D" k and oddness of 
7jfe w ^h respect to w to combine the two generators into 
this one expression. 

Many terms in this sum are suppressed by selection- 
ruled resonance conditions. In the second term, for ex- 
ample, most of the X ja pX kl s ~ e -«(^ Q ,3+^)t oscillate 
so quickly, compared to the weak coupling between the 
chain and the bath, that they average to zero unless there 
is a near-resonance u a p ~ w«5 7 . The second term is thus 
effectively a sum over nearly-resonant pairs of allowed 
transitions from states \0) and I7) into \a) and \S) re- 
spectively. So is the third term. 

Meanwhile, the first and fourth terms are effectively 
sums over allowed transitions from states \a) near E a 



Two general features of these coupled first order dif- 
ferential equations are decoherence and relaxation. De- 
coherence is caused by transitions from \a) and |<5), 



dp a 



— ^2 ( XjapXkp a Fj k — Xj aa X k SS^j k 



3,k 

Xk aa XjssT" k + XkspXjpsT^fy = —T a s (27) 







(here the second and third terms vanished because of our 
specific selection rules), and is exacerbated by transitions 
from nearby states. Setting 5 — a we see the relax- 
ation dynamics, in that the population p aa is flowing to 
and from the ppp at the selection-ruled transition rates 
2'Y^j k Xj a pX k p a Y^Paa and 2^2- k Xj a pX k p a T" k ppp. A 
stationary (and thermal) balance is eventually reached at 



PPp/Paa 



e -hu afj /k B T ( at ^ich pomt tue 



populations' effect on the off-diagonal coherences p aS , 



3,k,0 



pss 



a—S 



^S0~ 
jfcP/3/3 ' 

1=0 



~ . -r^0 a ~ \ 
r jk P00 + T 3k Pace), 



1=P 



6— a 



(28) 

also vanishes). The decay of the coherences allows the 
relaxation to proceed to a thermal equilibrium p? = 
e -H s /k B T J £ s provided the network of selection-ruled 
transitions does not isolate any subspace(s). For isolated 
networks each subspace will obtain its own a stationary 
thermal balance of populations constrained by the total 
available initial probability to be in that subspace. 

Remarkably, the decay of p a $ is offset, as it is on the 
central diagonal, by the nearly-resonant pairs of tran- 
sitions into I a) and \S) from any matrix elements pp 1 
within, and/or nearly within, the off-center diagonal that 
includes p a s. Indeed, a population- like coherence flow is 
established between p a $ and pp~ ( that is primarily into 
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the lower energy subspace when the temperature is low; 
quantum oscillations which are decaying in one subspace 
can in principle emerge in a lower energy subspace. 

V. NUMERICAL SIMULATIONS 

We simulate the effects of the matrix element equation 
in a variety of scenarios for two- and three-SQUID chains. 
In all cases we choose J and B to set the chain close-to or 
far-from resonant pairs of transitions (we can run at the 
critical point because of its exactly-resonant pairs). That 
way we can and do discard the oscillating coefficients 
from the equation, since their effect would average to zero 
anyway. The matrix element equation becomes a coupled 
first order differential equation with constant coefficients 
which we numerically integrate^. 

We consider large SQUIDs, R — 10/mi with / = 3pA, 
and small SQUIDs, R = lOnm with I = O.lpA. In both 
cases we imagine them to be spaced 4R apart in a solid 
crystal with mass density p = 5g/cm 3 and sound veloc- 
ity 5km/s. We set the Heisenberg splitting at 1.0GHz 
= 8J/h for a two-SQUID chain and 1.0GHz = 6J/h for 
a three-SQUID chain. We then choose SQUIDs with fre- 
quencies of 2B/h = 0.5, 1.0, or 1.5 GHz (below, at, or 
above the chain's critical point). We set the temperature 
to be 1/5 the SQUID frequencies, i.e. k B T/h = 0.1, 0.2, 
or 0.3 GHz. 

Although larger SQUIDs are possible, these param- 
eters keep the photon-induced decoherence rates^, in 
the absence of shielding, well below our phonon-induced 
rates. But still these SQUIDs' 4R separation is large 
enough to discard the j ^ k cross terms. At the critical 
point there are certain transition frequencies lo = <C 
whose cross term rates D^r° are identical to the 
j = k rates, giving a partially collective coupling to the 
bath. On the other hand, the small SQUIDs' 4i? sepa- 
ration is small enough to achieve a collective coupling to 
the bath for all to. 

In either case the matrix elements of the interaction op- 
erators are used to calculate all the constant coefficients 
in the matrix element equation and we then proceed with 
the numerical simulations for any initial state p(0). 

VI. DISCUSSION 

The relaxation and decoherence will thermalize the 
SQUIDs when the network of selection-rules is com- 
plete. For example, for two large SQUIDs dissipat- 
ing independently, the matrix elements of the interac- 
tion operators show that the network of selection-rules is 
\tpi) 1^2) «-> 1^4) and <-> IV's) <-> I'M- There is a 
reasonable pathway from any eigenstate to any other and 
these examples lead to thermalization. However, in the 
case of small SQUIDs close together there is a collective 
coupling J z and the network (for two SQUIDs) becomes 
l^i) \ip3) ^ 1^4) while 1^2) is a protected subspace. 



We attempt to summarize with a schematic: 

4 4 4 

3 ind 2 3 col 2 3 cri 2 

\ S \ \ 

1 1 1 

in which the first and second networks are those of in- 
dependent and collective dissipation. The third network 
is at the critical point lo\2 = for which the network is 
\i>i) iVte) 1^4) ^ l^)- Being at the critical point 
severs the \tp\) <-> \ifi2) link, as was done in the collective 
case, for uj\2Tjk = < 1, but not the \ip2) \fpi) link 
for which uJuTjk 3> 1. The only allowed transitions out 
of I -01 } and \4>2) are suppressed by the cold bath which 
is loath to supply the necessary energy. This helps to 
protect the population of ^2) and the coherence p\2- 

VII. CONCLUSION 

A chain of a few coupled SQUIDs exchanging their an- 
gular momenta with a phonon bath can be studied, in 
the Born-Markov approximation, with master equation 
techniques from the quantum Brownian motion model. 
The damping and diffusion give a matrix element equa- 
tion showing decoherence, relaxation, and the possibility 
for decaying quantum oscillations to emerge in a lower 
energy subspace. The relaxation adjusts the populations 
of the eigenstates and undermines their support for any 
superposition (coherence) between them, leading to de- 
coherence. 

The cascade of populations can occupy entangled 
states of intermediate energy, resulting in a surge of the 
entanglement of formation that indicates the number of 
singlets needed to form, from local operations and classi- 
cal communication, an ensemble of SQUID pairs p(t)r^. 
The entanglement is induced even though the SQUIDs 
are dissipating independently. 

The level spacings in the Heisenberg SQUID chain in- 
clude pairs of resonant transitions which are necessary 
for coherent oscillations to decay into a lower energy 
subspace where they can decohere more slowly. In this 
phenomenon, a superposition of two eigenstates relaxes 
coherently into a superposition of two lower-energy eigen- 
states with the same energy difference as the upper two. 

Small SQUIDs close together exhibit a collective cou- 
pling to the bath which can give a protected subspace 
and enhanced or suppressed transition and decoherence 
rates. In effect, the network of selection-ruled transitions 
is broken into isolated pieces. When the level spacings 
cooperate to allow coherence flow in these sufficiently iso- 
lated pieces, decaying quantum information can reappear 
and be sustained in a lower energy subspace; this is the 
idea behind "automatic quantum error correction"—. 

Another feature of the Heisenberg SQUID chain is the 
critical point level crossings where an allowed transition 
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vanishes along with its frequency. The network of 
selection-rules in effect acquires some features of the 
collective behavior as uoTjk — > <C 1; even large SQUIDs 
spaced well apart, when tuned to the critical point, 
can have extended coherence times. More generally, 
this suggests some qubit encoding schemes might be 
augmented with intcr-qubit couplings to obtain some or 
all of the benefits of a collective coupling to the bath. 
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APPENDIX A: COHERENCES AND 
POPULATIONS 



the Fourier transform of the <fr within the ring, we have 

ks 

In the thin ring approximation, \k\Vb <C 1, the Fourier 
transform becomes 

<% -j2irRbJ 1 (k x R)n£, 
where J\(k x R) — Ji(|fc|.Rsin#) is the first order Besscl 
function. The polar angle 8 is the angle between the z- 
axis of the ring and the wavevector k while _L k lies in 
the plane of the ring. Choosing ei(fc) to lie in the plane 
of the ring, i.e. ei(fc) = n^, we obtain 



The instantaneous state of a spin-chain can be de- 
scribed by a density operator p(t) = J2 n p n \ip n ) {ip n \ 
which averages "outer products" of pure states \ip n ) (ip n \, 
weighted by their probabilities p n , into a statistical 
mixture. In the energy eigenbasis, a diagonal element 
(a\ p\a) is the probability of obtaining the eigenstate \a), 
sometimes called the population of \a). An off-diagonal 
element (a\ p \(3) results from including superpositions of 
eigenstates, e.g. \ipi) = c a \a) + cp \(3); it has an evolv- 
ing phase, (a\p(t)\0) = e i ( E «- E ^ t / h (a\ p{0)\f3), which 
indicates the chain's coherent dynamics. For this reason, 
these off-diagonal terms are called coherences. 

APPENDIX B: THE CHAIN-BATH COUPLING 

The first SQUID has a current Iaf confined to its ring 
of cross-sectional area b: j\ — (7 '/b)af<j> within the ring, 
ji = elsewhere. Here <j> is the azimuthal unit vector in 
cylindrical coordinates centered on the ring. With 

Vi = (m e /e) J d 3 rj! ■ u 

and 

p2w pR+Vb/2 rVb/2 

d(j) drr dzcf)e~ lk ' i ' 

JO JR-Vb/2 J-Vb/2 

I 



V 1 = \a{ Y,\f\k 7 \RJi(k x R) 4i) ' 

k 

with coupling constant A = 27r/(m e /e) Rhc± /2pV. 

The other SQUID rings are centered not at f = but 
are evenly spaced, a distance d apart, along the x-axis. 
The analysis for each SQUID's coupling to phonons is 
calculated in its own coordinates fj centered at Xj = 
d(j — l)x so that f=Xj + fj. The creation and annihila- 
tion operators' phase factors e ±lk ' r become e ±lk ' Xj e ±lk ' r i 
while the rest of the calculations, in the fj coordinates, 
are exactly the same as before. The total coupling is thus 

3 k 



APPENDIX C: THE BATH CORRELATOR 

To evaluate the bath correlator we perform the trace 
in the bath's energy eigenbasis by summing over the di- 
agonal matrix elements: Tr B [(f>j(t) <j)k(t') Pb] = 



]T \J\k\R\k'\R.h(k x R) Ji(k' x R) (n| (a % SV tS > + ^K^') (a^ + at^'K^) & ~ " - \n) 

n k.k' 



\k\Rjf (k^R) 



f°° dujuj 2 f 2n ' * * 

(U/8tt 3 ) / —g— / d8 sin 9 dcj>ujT R Jl{LOT R sin 9) 

JO C J_ Jo Jo 



coth( 



tlLO 

2k B T 



cos(wr) — i sin(wr) 



cos(uiTjk sin 9 cos <f) . 



(JV„ + l)e- i "' T +iV w e i( 



±ifc-(x 
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The steps leading to the second line are as follows. The 
double sum Ylkk' c °U a P ses to a single sum (which 
we convert to a fc-space integral in spherical coordinates) 
because the only non-zero cross-terms ~ fifcfcf • They are 



(n| a %x e 
and 



a bi e 



-ik' -x k 



-iujT ^ik-Xjk 



(n au e 
x 1 fei 



Hihk'^e 



where the factors e ±luJT arise from being in the inter- 
action picture. Here we have abbreviated u\{k) = uo 
for the angular frequency of the mode with wavevec- 
tor k and transverse-in-plane polarization and written 
Xjk = Xj — Xk for the vector connecting the (j, fc)-pair 
of SQUIDs. We have also used tr = R/c±, which is 
half the time it takes a phonon to traverse a SQUID. We 
switch the order of the sums, J2 n J2k ^ SfeSn' an< ^ 
write = J2 n n k\ ~ — g~ ~ f° r the thermal-average oc- 
cupation number. It sums to N u = i/( e W=BT _ 

E -hujn Ul /k B T tt ^ -huJini/k B T 

E e -fia,n Sl /fc B r n E e-^n./fcBT " 

kl . ' . 
v v ' 

l/(e» u /*fl T -l) 1 

Now every function in the bath correlator besides the 
e ±ik-x jk - 1S an even function of k. The sum over wavevec- 
tors thus selects the cos(fc • xjk) part of e jk . Be- 
cause the (j, fc)-pair of SQUIDs are positioned on the 
x-axis, we use k x = \k\ sin 9 cos <fr to obtain cos(fc • Xjk) — 
cos(uTjk sin 9 cos</>), with Tjk = d(j — k)/cj_ the phonon 
transit time between the SQUIDs. Finally, it can be 
shown that {2N 0J + 1) = coth(fkv/ '2k bT) so that we can 
write 



N„ = 



{N u + l)e 



=coth( 



tlLO 

2k B T 



) cos(wr) — i sin(wT). 



APPENDIX D: THE SPECTRAL DENSITIES 

We can write the bath correlator as an integral over 
bath frequencies, Tr B[<j>j(t) <pk(t') Pb] = 

f°° Tiuo f°° 

I du> Jjk(u) coth(— — — ) cos(ujt) — i / du> Jjk{u) sin(o;r), 
Jo 2k B T J 

v ' S v ' 

characterized by spectral densities 

for which we need the angular integration 

<djk(oJTji) = / d9 sin 9J 2 (lotr sin 9) 2wJo(u)Tjk sin 9) . 

Jo . ' 

Jd<p cos(oJTjfc sin 8 cos 0) 



For the single-SQUID (j = k) 9 integration we obtain 



J. 



n 



, . hm 2 I 2 R 4 Lu 5 , r 3, r 5 , , n9n 



independent of j and with p F q the generalized hyper- 
geometric function. Although we use these exact " p F q " 
spectral densities in our numerical simulations, it is help- 
ful to know that for a "small ring" {lotr <C 1) the single 
small SQUID spectral density is approximately 



2\2 



2Hm 2 e I 2 (irR 2 ) 
5A 2 7r 2 e 2 4pc^ 



while for a "large ring" (wtjj 3> 1) the single large SQUID 
spectral density is approximately 



2hm 2 I 2 R 
( 

A 2 e 2 Apc 2 ± 



Next, for uiTjk < 1, as can happen for "large" SQUIDs 
at a level crossing where ui w or, regardless of ui, for a 
few small SQUIDs spaced only a few small SQUID radii 
apart, the cross term spectral densities must be consid- 
ered. It turns out that to first order in uiTjk we can use 
the single SQUID j = k rates for the j ^ k cross terms, as 
follows, lot jk <C 1 implies utjj <C 1 so we reconsider our 
9 integrations with J\ {ujtr sin 9) w ^LJTRsmO to obtain 



2-7T J d9 sin 3 Jo(ujTjk sin t 



Ojk(wT R ) 

/4 



,2^2 



47T 



UiTjk cos(u>Tjk)- (l-w 2 T 2 fc ) Sm(u>Tjk) 



UJ 3 T 3 , 



which, to first order in cuTjk, give spectral densities iden- 
tical to those of a single small SQUID 



Jj% J)T^) when ^ fe « 1 

and we may as well use the equivalent p F q spectral den- 
sities. 

Finally, we assume that when uiTjk 3> 1 we can neglect 
the cross term (j ^ k) spectral densities. The assump- 
tion relies on the Jo(oJTjk sin 9) kernel of the 9 integration 
oscillating quickly enough between positive and negative 
values that the integral never accumulates any significant 
value. 

Again, the p F q spectral densities are the ones we use 
for j = k and for close-together or far-apart SQUIDs 
satisfying uiTjk <C 1. For far-apart SQUIDs satisfying 
ujTjk >lwe use Jjjtk(w) = 0. 



AS) 



APPENDIX E: CONSTANT (MARKOV) 
COEFFICIENTS 

The coefficients of diffusion, rcnormalization, anoma- 
lous diffusion, and damping are obtained by convolv- 
ing a coswr or sinwr with the Vjk(T~) or /Xjfc(r) ker- 
nels. But the only time-dependent terms in the inte- 
grands for the kernels are cos lo't or sino/r. Since these 
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time-convolutions are relatively simple, in our calculation 
of the coefficients we switch the order of the integrals, 
J dr J dui' =4> J duj' J dr. 

For the diffusion and damping we find that the inte- 
grations Jg dr cos uj't cos uit or J Q * dr sin uj't sin wt are 



sin(o/ — uj)t sin(o/ + u>)t 



which behave, as t — ► oo, like Dirac delta functions: 
^ [£(u/ — w) ± + u>)]. This is because they oscillate 
with lo' at a frequency t everywhere except at ±w, where 
they spike to a height ~ t and width ~ 1/t. As long 
as the frequency of these oscillations is much faster than 
any features of the spectral densities (i.e. t 3> t^tr), 
of the cutoff (i.e. t > 1/A), and of the hyperbolic cotan- 
gent (i.e. t » Ti/ksT) then the only contribution to the 
ui' integral comes from the spikes at ±u). Subject to these 
rough criteria, 

tyz>T jk ,T R t>l/A t^>fr/k B T, 

the coefficients of diffusion and damping approach con- 
stant values which, thanks to the Dirac delta functions, 
are easy to identify: 



and 



1%= ~Jjk(u)= lim<yft(i). 



For the renormalization and anomalous diffusion, the 
integrations f* dr sin uj't cos ujt or J Q * dr cos o/r sin o>t 
are 



n 2 a'+m 

2 .± 



,2 ui' —ui j 



Uj' + W 



whose behavior for t — > oo is not so clear. However, for 
large t these terms oscillate so fast with lj' that they 
too average away the time dependence and we can at 
least define, using an explicit exponential cutoff / A , 
constant values for the renormalization, 



/•OO 

r%= lim / doj'J jk (J)e-"'l A 



2 (jj'+uj 



t sin 



4 



w' + ui 



and anomalous diffusion, 



A 



jk 



dui'Jj k (uj') coth( 



2k B T 



rt— >00 

/ A / c?t cos w't sin ur, 
Jo 



which we use symbolically without actually ever evaluat- 
ing them. 
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